Reduced coherence in double-slit diffraction of neutrons 
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In diffraction experiments with particle beams, several effects lead to a fringe visibility reduction 
of the interference pattern. We theoretically describe the intensity one can measure in a double- 
slit setup and compare the results with the experimental data obtained with cold neutrons. Our 
conclusion is that for cold neutrons the fringe visibility reduction is due not to decoherence, but to 
initial incoherence. 
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I. INTRODUCTION 

We provide a theoretical description of the inten- 
sity pattern in double-slit experiments with neutrons, 
with specific attention to the cold neutron diffraction 
(A w 20 A) carried out by Zeilinger et at in 1988 Q. 
The result we obtain is shown in Fig. [TJ 

Usually, the main problem in the analysis of diffraction 
experiments is to establish exactly which causes bring 
about the reduced coherence experimentally inferred 
from the detected signal. This effect can be produced 
both by the initial preparation of the beam (namely, the 
non-dynamical incoherence) and by the interaction of its 
constituting particles with the environment (namely, the 
dynamical decoherence) @, |, i, H [1, 0, B S, El EI • 

We find that in the experiment of Ref . [l[ decoherence 
does not play any role in the fringe visibility reduction, 
which indeed is entirely due to incoherence of the source. 
This conclusion is opposite to that of Ref. [ll| , where it 
is claimed that decoherence is essential for explaining the 
data from 

We provide calculations and numerical simulations in 
support of an unexpected incoherence cause, which we 
propose plays a role in explaining the experimental data 
of Ref. [1]. This cause is that the width of the wave 
function impinging on the double-slit is comparable with 
the size of the double-slit setup (i.e., the distance be- 
tween the two slits and the slit apertures). We argue 
that this feature of the incoming wave function leads to 
a slight difference in transverse momentum between the 
two wave packets emanating from the grating, such that 
the centers of the packets move apart. This momentum 
difference was already suggested in Ref. [l]J , but no phys- 
ical explanation was given. It proves relevant for fitting 
the data of Ref. Even though it does not reduce 
the fringe visibility, it can significantly affect the shape 
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Figure 1: Comparison between theoretical prediction derived 
here (solid line) and experimental data taken from [l[ (solid 
circles) for neutron double-slit diffraction. 



of the interference pattern, changing the overlap of the 
two packets and so the position of secondary minima and 
maxima. 



II. 



VISIBILITY REDUCTION IN THE 
INTENSITY PATTERN 



We now describe the various causes of reduced coher- 
ence which can occur within neutron interferometry. Let 
the y direction be the direction of propagation of the 
beam. We assume that the grating is translation invari- 
ant in the z direction, and that the motion in the y direc- 
tion is essentially classical, so that the problem reduces 
to one dimension, corresponding to the x axis. 



Initial incoherence 
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Concerning the non-monochromaticity, every spectral 
component of the beam contributes incoherently to every 
other. In our case, assuming that the detection screen is 
parallel to the x axis, the intensity observed on the screen 
at coordinate x is 



I(x) = / d\f(X)h(x), 



(1) 



where /(A)dA is the wavelength distribution of the beam 
and I\(x) the intensity corresponding to a single wave- 
length A. Thus, one can continue the analysis with a 
single wavelength and postpone the integration ([lj to 
the final step (see Sec. IIIip. 

Concerning the collimation of the beam, a relevant 
cause of incoherence is that the particle source emits ran- 
dom wave functions, whose randomness lies in the devia- 
tion of its direction of propagation from the y direction, 
corresponding to imperfect collimation. This cause of in- 
coherence will be taken into account in our model of the 
beam by means of a random wave vector k. 



B. Decoherence 

In the treatment of decoherence, we make explicit use 
of the fact that we are concerned with neutrons, in par- 
ticular at the low energies of Ref. [l| . In this case, in fact, 
the only relevant decoherence channel — if any — consists 
in collisions with air molecules, so that the dynamics can 
be modelled within a Markovian description of the scat- 
tering event, in particular in the large scale approxima- 
tion allowed by air molecules (see @, H, 03] ) • The model 
obtained in this way allows an estimate of the neutron 
coherence time r co h [H, [3] : 



1 



7"coh 



(2) 



where P(Qg) is the environmental pressure at the tem- 
perature Of, CTtot the total cross section of the scat- 
tering events, fcs the Boltzmann constant and m a ; r ~ 
4.8 • 10~ 26 Kg the mean mass of air molecules. (Eq. @ 
takes into account a correction by a factor 2ir, usually 
missing in the literature, that was theoretically predicted 
in [H,Il3] and experimentally checked in [§].) 

The result is that r co h is much greater than the time of 
flight, even in extreme experimental situations. In fact, 
even considering a surrounding pressure of 1 atm (though 
typically "the beam paths along the optical bench" are 
"evacuated in order to minimize absorption and scatter- 
ing" [l|) and room temperature, for an estimated total 
cross section of 10~ 27 m 2 [lj| we obtain that T con ~ 140 s. 
This time not much smaller than the neutron lifetime — 
shows that the coherence is fully kept for the duration of 
most experiments. For instance, in [l[ the mean time of 
flight is T ps 0.023 s, several orders of magnitude smaller 
than T coh . 

For this reason, in the following we shall regard the 
neutron beam as uncoupled from its environment and 



shall describe the dynamics through the usual unitary 
Schrodinger evolution. 

This conclusion contradicts that of Sanz et al. in 
Ref. [IH , who claim that "decoherence is likely to exist in 
Zeilinger et al.'s experiment." Indeed, their own calcula- 
tions do not support their conclusion. The basis of their 
claim is that the damping term A t in their expression for 
the observed intensity pattern [see their Eq. (27)] turns 
out, when fitted to the data, to be nonzero. However, 
this term could just as well represent incoherence instead 
of decoherence (as does for example the similar quantity 
A in their Eq. (16), or the coherence length in Eq. (47) 
of Ref. [1]) [2(J. Indeed, our estimate of the coherence 
time shows that the damping A t cannot be attributed to 
decoherence. 



III. MODEL FOR THE INTENSITY PATTERN 

In order to obtain a theoretical description of the in- 
tensity pattern measured in the experiment, we precisely 
describe the whole evolution of the neutron beam, from 
its production to the arrival on the detection screen. 

We set up a concrete model for what the wave functions 
of the neutrons in the beam look like, and thus obtain 
the densit y m atrix and the intensity pattern. (Another 
approach [131 ] is to guess the density matrix from infor- 
mation theoretic principles, but we prefer to avoid the 
invocation of such principles.) 

The width A of the entrance slit fixes the spatial ex- 
tension of the (random) wave packets. The wave packets 
tpo produced by the source are modeled, when passing 
the entrance slit, as Gaussian wave packets with mean 
and standard deviation sq — A/^/l2, i.e., the standard 
deviation of the uniform probability distribution over the 
interval of length A. Thus, 



ip (x 



x 2 

(x) ~ exp I ~ + ikx ) 

V 2sq / 



where we take the random wave number k to have a 
Gaussian distribution with mean and standard devi- 
ation tTfe. As established in Sec. Ill Bl away from the grat- 
ing, the neutrons follow the free Schrodinger evolution, 
so that — after the flight over the distance Lq — the packet 
just before the grating (the incoming wave function) has 
the Gaussian form 



ipin(x) ~ exp I - (1 - fy) 



8(k)] : 



2s 2 



+ ikx 



with 7 = AL /(47tsq), s = s y/l + 7 2 and S(k) = 
LoXk/ (2ir). In Sec. lIVI we shall present two different mod- 
els for the wave function immediately after the grating 
(the outgoing wave function ip ou t) a s a function of the 
incoming wave function tpi n . 

We now consider the free propagation of the outgo- 
ing state ijj out toward the screen. We obtain the density 
matrix p out by averaging |V'out)(V'out| with respect to the 
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probability distribution of the wave number k. The in- 
tensity measured on the distant screen is proportional to 
the diagonal elements of the density matrix pr at the 
time of arrival T — mL\/(2irh) (with m the mass of the 
neutron) — as already stated and well motivated in the 
literature [1, [13]: I\ (x) ~ pr(x,x). Before identifying it 
with the measured intensity, we have also to consider the 
integration ([I} over the wavelength and the finite spatial 
resolution xq of the detector. For the latter, we shall 
consider flat response on the interval [x—Xq/2 , x+xq/2] 
around each position x, so that the intensity is finally 
given by 

i r px+xo/2 

I(x) = — d\f(X) dy I x (y), 

x J Jx-x /2 

where xq and /(A) are directly adoptable from Ref. 



IV. PASSAGE THROUGH THE GRATING 
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Figure 2: Comparison between the theoretical prediction de- 
rived with our first model for the passage through the grating 
(solid line) and experimental data taken from [1[ (solid cir- 
cles). 



We describe two models for the passage of the neu- 
tron through the grating. The first is a simplified model, 
which can be treated analytically and features the mo- 
mentum difference between the two transmitted packets 
in the transverse x direction. The second, on which Fig. [I] 
is based, is more flexible and achieves better agreement 
with experimental data. 



A. First model 

The central assumption of this model is that the out- 
going wave function is of the form 

lpout(x) = F(x)lp in (x), 

where F is a complex function whose modulus is < 1 
and signifies the degree of transmission, while the phase 
signifies the phase shift during the passage through the 
grating. The simplest model assumption would be to set 
F(x) — 1 for x € [x\ — ai/2, x\ + a\/2] U [x? — »2/2, X2 + 
02/2] and F(x) = otherwise, where the two intervals 
represent the two slits, their apertures a%, ai and distance 
d are known and Xj = (— l) J (aj + d)/2 for j = 1,2 
is the center of the j-th slit. With this choice of F, the 
multiplication by F is a projection operator. However, 
to enable analytical computations we find it more useful 
to set 



F(x) 



1 



■ exp 



(x - £i) 2 



2(7? 



1 

cr 2 



■ exp 



x 2? 



2a\ 



where Oj — aj/y\2. 

Now I\(x) can be calculated exactly. In Fig. [2] the 
comparison between the theoretical intensity which fol- 
lows from this model and the experimental data of Ref. [l| 
is shown. Apart from a constant background subtracted 



from experimental data, there is only one fit parame- 
ter <7fe, for which we have found o~k = 10 5 m/s. This 
value corresponds to a x— velocity with standard devia- 
tion equal to 0.006 m/s, while the mean y— velocity dis- 
tribution is 214 m/s [l|. 

It is useful to look at the explicit formula for the out- 
going wave function. For simplicity, set k = 0. Then 



Ipout(x) = 1pout,l(x) + l/W,2(z), 



where 



iPjZ 



4>aut,j( x ) ~ ex P [ - ( X ~ &) ( a 3 - 2^2) + 

and where we have introduced the following quantities, 
for j = 1,2: a, = l/(2s 2 ) + l/(2a 2 ), = x j /(2a j ^) 
and pj = h'jXj/(a'j + s 2 ). Therefore Vwt is the sum 
of two Gaussian packets with momentum in the x di- 
rection equal to p\ and p 2 - Moreover, since x\ and xi 
have opposite signs, the same holds for £1 and £2 and, 
consequently, for p\ and p2- In other words, the pack- 
ets produced by the grating have two opposite momenta 
in the transverse direction, as a consequence of the finite 
width of the neutron wave function compared to the scale 
fixed by the grating size. (Note that if tp- ln is very spread 
out, i.e., if s is large, then pj ~ (ctJ + s 2 )" 1 is small). 

Also in the presence of a common drift expressed by k, 
the packets produced by the grating have wave number 
in the transverse direction equal to k+pi/h and k+p2/h, 
with analytical expressions for the wave functions and pj 
more complicated than before. 

Even though the agreement with experimental data is 
not completely satisfactory, this model illustrates the ori- 
gin of the transverse momenta p\ and p 2 that, as already 
noticed by Sanz et al. in Ref. [Tl| . move outwards the 
x— position of secondary minima. 

The values pj can also be predicted numerically by 
solving the Schrodinger equation for the passage of the 
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neutron through the grating, treating the problem as two- 
dimensional with the double-slit modelled as an impene- 
trable potential barrier. Our numerical simulations were 
not accurate enough to yield realistic values, but did con- 
firm that the pj are effectively nonzero, point in opposite 
directions (outwards), and have roughly equal modulus. 
Numerical simulations also suggest that the values of pj 
are substantially the same for all A within the range of 
wave lengths selected for the experiment. 

B. Second model 

The second model allows to fit p\ and pi to the data. 
It is based on the assumption that the outgoing wave 
function is of the form 

i>out(x) = ci <fi{x) + c 2 <pi{x) , 

where Cj (pj is the wave packet emanating from the left 
(j = 1) or right (j = 2) slit, which we assume to be of 
Gaussian form: 

r„./Vl — - -(x- Xl ) 2 /4a^ ikjx 

ft(l, "(2^)i/ 



allowing for the possibility that the drift kj and the 
weight Cj are not equal for the two packets, where <jj and 
Xj are the same quantities as introduced in Sec. IIV Al 

We write the wave numbers kj as kj = k + Pj/fr, i-e., 
the sum of the wave number k of the incoming packet and 
the transverse momentum pj acquired during the passage 
through the grating, as previously highlighted. 

We estimate the coefficients Cj as 

Cj = Cj {k) = 1^(^)1 . 

Thus the wave packet emanating from slit 1 has greater 
weight | ci | 2 if the incoming wave packet is closer to slit 
1. (For further discussion see Ref. [Hj].) 

Fig. Q] shows the resulting prediction of this model. 
The following parameters were obtained by means of a 
fit procedure: the two velocities px/m — — 0.0034 m/s 
and pijin = 0.0029 m/s 21], and the angular divergence 
of the beam ak = 4976 m -1 (which corresponds to a 
x— velocity distribution with standard deviation equal to 
0.0003 m/s). Moreover, also in this case a constant back- 
ground has been subtracted from experimental data. The 
Pj have been assumed to be the same for all wave lenghts 
A selected for the experiment. 
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To our knowledge, direct values of otot under the consid- 
ered conditions for the molecules which compose the air 
are not reported in the literature. Nevertheless, starting 
from related measurements [i"7l ] and using typical trends 
of the cross section, it is possible to estimate very roughly 
the bound <7tot 10~ 27 m 2 . However, even if a correct es- 
timation of otot led to a greater value, r co h would still be 
much larger than the time of flight T since five orders 
of magnitude lie between our estimate for r co h and T. 
Furthermore, we have assumed very unfavorable experi- 
mental circumstances. 

In particular, the caption of Fig. 2 in ,11] is incorrect 
in so far as Fig. 2(a) does not include incoherence, and 
Fig. 2(b) may represent either incoherence or decoher- 
ence. 

These values of pi and P2 agree qualitatively with the 
estimates in Ref. [Til ]. Nonetheless, the estimation in [ll| 
is dubious, as it is based on an incomprehensible appli- 
cation of the Heisenberg uncertainty principle. 



